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ABSTRACT 

The turbulent boundary layer under a freestreani velocity that 
varies sinusoidally in time around a zero mean is considered. The 
flow has a rich variety of behaviors including strong pressure gra- 
dients, inflection points, and reversal. A theory for the velocity- 
and stress profiles at high Reynolds number is formulated. Well- 
resolved direct Navier- Stokes simulations are conducted over a 
narrow range of Reynolds numbers and the results are compared 
with the theoretical predictions. The flow is also computed over 
a wider range of Reynolds numbers using a new algebraic turbu- 
lence model; the results are compared with the direct simulations 
and the theory. 

NOMENCLATURE 

f\y fi unknown nondimensional functions, see Eq. (3), (4)... 
L == r 0 /u> outer length scale 
Re - l 'o Si / v Reynolds number 
t lime 

u, t\ u‘ velocity components 
u* - ti V fluctuating velocity component 
C - < u > mean velocity component 
Ut laminar mean velocity component, see Eq. (3) 

Uo peak freestream velocity, see Eq. (I) 
l r ao freestream velocity, see Eq. (1) 
u* peak friction velocity, see Eq. (5) 
u r friction velocity (can be negative) 

U* = U / u T wall units 

z, y, z st Teamwise, normal, span wise coordinates 
y + El y\u r \iv wall units 

S\ = y/2v/w laminar boundary-layer thickness, see Eq. (3) 

6 turbulent boundary-layer thickness, see Eq. (5) 

* K arm an constant 

A,, A r periods in z and z directions 

v kinematic viscosity 

<t> = utt phase angle 

<fr% phase shift, a function of Re y see Eq. (7) 
w frequency, see Eq. ( 1 ) 
r = U y — < uV > total shear stress 
< > average in z, z and/or ensemble 

INTRODUCTION 

The flow under consideration is the turbulent version of a classi- 
cal exact solution of the Navier-Stokes equations, known as Stokes’ 
second problem (1). The domain is the half-space over a station- 
ary flat plate and the freestream velocity is given by 

Uoo(t) = l r Q cosM). (1) 


The three parameters (peak velocrty f/* 0 , frequency and kine- 
matic viscosity v) combine into a single nondimensional parame- 
ter; a Reynolds number Re. The laminar solution is 

Ui(y,t) = tfco(<) - U t e->' 6 ‘ cos M - y/6 t ) (2) 

where 61 is the laminar boundary-layer thickness. 

This flow is attractive because of its rich time- dependent be- 
havior, and because it is naturally homogeneous in the directions 
parallel to the plate (this is true only because the mean value of 
Uooit) is zero). This justifies imposing periodic conditions in a di- 
rect simulation and considerably improves the statistical sample. 
In addition, the Reynolds-averaged quantities are functions only 
of the normal coordinate and the time. This provides a good test 
case for turbulence models, in which their ability to treat unsteady 
and inflectional velocity profiles can be tested with little numerical 
effort and with high numerical accuracy. It also exhibits reversal 
of the velocity without the usual viscous-inviscid coupling diffi- 
culties and singularities encountered in spatially-separating flows. 
The fully-developed flow will most likely be close to the common 
turbulent boundary layer, except that the lack of entrainment will 
prevent the formation of a sharp laminar-turbulent interface. 

The known results concern primarily the transition Reynolds 
number, Rc ss 600, and the stability of the laminar solution 
to small disturbances (2). Detailed experimental results are not 
available for validation. However the only source of error is nu- 
merical: an approximate treatment of the boundary- layer growth 
is not needed here as it was in Ref. 3- The numerical errors can 
be reliably estimated by conducting simulations with improved 
resolution, enlarged domains, and so on, and by using the expe- 
rience gained while simulating other flows (2, 4). The agreement 
with well-established laws (e. g., the log law) also helps build 
confidence in the numerical results. 

TENTATIVE THEORY OF THE FLOW 

The theory is tentative in that it is based on assumptions which 
are plausible and consistent with current thinking, but cannot be 
rigorously proven. As is often the case in turbulence, only neces- 
sary conditions will be obtained. Also, the evidence supporting 
the theory is encouraging, but still limited. We focus on the mean 
velocity U[y,t] and the total shear stress r[y,f]. They are linked 
by the momentum equation, U t - ~ t v . For clarity we use the 

following notation: square brackets surround the arguments of a 
function, while parentheses group terms in an expression. Most 
of the derivation can be generalized to other laws than cosiu/tj for 

Uoo/Uo. 
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Define the phase <t> and the length scale L. Dimensional analysis 
shows that U and r must have the following form: 


U[y,t\ -M‘) = 0./,[J,* t (3*) 

r[»,<] = Uih[*,4,Ik] (3b) 

for some nondimensionaJ functions /j and / 2 . Recall that the 
averages are taken over x and z with sufficiently large domains 
and by phase- averaging. The momentum equation and boundary 
condition become 


dh = dh 

d<t> d{y/L )’ 


/ 1 [0,^Ae] = (4) 


So far no approximations have been made. The functions f\ and 
/ 2 are functions of three variables; the objective is to separate the 
dependence on Re from the other two as much as possible. 

Outer reg ion 

The Reynolds- number-similarity hypothesis (5) asserts that, 
away from the wall, the gross features of the turbulence (at least 
U and r, which do not depend strongly on the small-scale motion) 
can be independent of the Reynolds number provided that proper 
velocity- and length scales u m and 6 (different from Uq and L) are 
chosen to normalize them. In the outer region we replace Eq. (3) 
by 


i%,<) -(/«[<] = u- /,[!>.&]. 

(5a) 

r[y,<] = u ,J / 4 [|,*,*e]. 

(5b) 


This equation is just as general as Eq. (3), since u* /U q and 6/ L 
are functions of Re, but / 3 and / 4 are now of order 1, as is the 
range of y/6 over which most of their variation occurs. We can 
define u m as s /r TTiaz where r maz is the peak value of the shear 
stress at the wall during the cycle. For 6 there is no such obvious, 
quantitative definition (especially since the flow does not have a 
sharp laminar- turbulent interface). Since the time scale is 1/w. 6 
must be proportional to «*/«;, We can use this to define 6. We are 
free to choose the value of the ratio, and we choose 1. Therefore: 


Clearly when the wall shear stress reverses (r(p,t] vanishes). Hus 
assumption breaks down. The theory will n<>i apply near the phase 
of reversal, but this may not have much effect on the rest of the 
cycle. 


In the overlap region both Eq. 
to be valid and we obtain 


(7a) and Eq. (8) are assumed 


? few * 


y |u r | u m 6] 


To exploit the overlap argument we can, for instance, vary Re in 
Eq. (9) while holding y/6 and 4> + <f> 0 fixed. We obtain 

Uoo d (Uo\ Uo d /tfooX d<fro = u r df 7 rf(log[tt**/i/j) 

Uq dRe\u m ) u* d<j> V U 0 ) dRe u* d(logfy+]) dRe 

(10) 

The only term that depends on y is the last one, which contains 
d/r/d(log[y ■**]). Therefore this quantity must be a constant, e.g., 
1 /*; again the overlap argument predicts a log layer 

(/t = log[y + ]/ k + C). 

The reason d> 0 (R«) was introduced is not apparent. Consider 
the variation of the various terms in Eq. (10) as a function of 4>. 
If <f> o were independent of Re , the second term would be 0. Then 
during the cycle, u r /u* (the third term) would be proportional to 
Uoo/Uq (first term), so that the wall shear stress would have to 
have the same phase as the freestream velocity, which is physically 
not correct. The variable-^ model allows the deceleration to enter 
through d(U qq / Uq) / d</> (second term) and make the wall stress 
“lead” the freestream velocity and reverse earlier (if d4> 0 /dRe < 0). 

Dependence on 6 and Sn 

We conclude the analysis by exploiting the dependence of Eq. 
(10) on 4> and <po- Using Eq. (1) we see that the left-hand-side is 
proportional to cos(^ + 4>\) with <t>\ defined by 

tan($i] = d<t>o/d(\o%\u* /Uo}). The right-hand-side is proportional 
to u T /«' which is a function of <t> 4- ^o- Therefore </>\ - 0 O is a 
constant. Since <t > 0 is meaningful only up to an additive constant 
(see Eq. (7)), we can choose the constant so that <j> 0 ss 0,. Con- 
sequently tan[<£ 0 j = dd>o/d(\og[u m /U 0 \), which integrates to 


*=-. (6) 

u? 

Thus 6 is proportional to the thickness of the turbulent region, but 
not equal to it. Because of Eq. (6), fz and / 4 satisfy an equation 
analogous to Eq. (4(a)), simply replacing L with 6. 

We now introduce the major assumption: the only dependence 
of f z and / 4 on Re is a phase shift 4> 0 [Rej (as we shall see, the 
theory is too limited if we assume that / 3 and / 4 do not depend 
on Re at all). The idea behind this assumption is that at different 
Reynolds numbers U — f/* and T (scaled with u* and 6) may have 
the same behavior, but not at exactly the same phase relative to 
Uqo. The assumption is expressed as follows: 

/j = /s[|*<A + ^°[Re]], (7a) 

ft = + ^o!^ e lj- (7b) 

Note that / 5 and /* are functions of only two variables. On the 
other hand we have introduced three new functions of Re: u m /Uo, 
6/ L and <p Q (but recall that the first two are related by Eq. (6)). 

Inner region 

Here we apply the usual arguments. The variation of r with y is 
neglected, so r[y,<] ^ r[0,<], and the momentum equation is not 
considered. Define u T = ^/r[0,(] where the square root is taken 
to have the sign of the argument. The velocity U is assumed to 
follow the law of the wall: 

V = (8) 

L v J 


sin{d>o] = 



( 11 ) 


where A is a constant. Also, u* was defined so that the peak value 
of u r /it* is 1. Therefore 


^=cot[* + *o]. (12) 

Finally, rearranging Eq. (10) and using Eq. (6) to eliminate 6 , 


Ak 

2 sin[<£o] 


+ cos 



d(j ) 0 

sinf^o] 


dRe 

Re 


(13) 


This integrates to 


k Uq . J , 

- — COS(0»] + log 



= log[Ae] - B 


(14) 


where B is a constant. It turns out that 4>q is small, of the order 
of 0.15 radians, so that cos[^ 0 ] 0.99. If it is replaced by 1, Eq. 

(14) involves only u*/Uq and Re, 

Equations (11), (12), and (14) provide strong predictions that 
can be checked against sufficiently accurate simulations. Such 
a check will also provide the values of A and B. The theory 
predicts that as Re — * oo, u m /Uo y 6/L , and <p 0 all tend to 0, 
although slowly (roughly like 1 /log(Re), as in other flows). Also 
note that the nondimensional shear-stress gradient dr+ /dy* is of 
order (Uq/u *) 2 / Re 2 \ it also tends to zero as Re — ► oo so that the 
“almost-constant-stress” assumption for the inner region becomes 
more and more valid. 
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The known limitations of the theory are that it applies only 
at sufficiently high Reynolds number ami that it does not apply 
near the phase of reversal of the wall shear stress. This failure near 
reversal is obvious since Eq. (12) predicts that r has a double zero 
when it reverses, which is counterintuitive and is not indicated by 
the simulation results. It may be possible to extend the theory to 
cover the reversal region by matching asymptotic expansions, but 
this has not yet been achieved. 

DIRECT-SIMULATION RESULTS 

The numerical method used to solve the time-dependent, three- 
dimensional, incompressible Navier-Stokes equations is described 
in detail in Ref. £. It is fully spectral in space and second-order- 
accurate in time. The method was further tested by solving the 
linear stability equations for the laminar solution using Floquet 
theory (in that case the time integration used a fourth-order- 
accurate, implicit, reversible scheme). The results obtained by 
Hall (7) with a very different method were reproduced, which 
shows that the spectral basis functions that are used are well 
adapted to the flow. The results indicate linear stability at all of 
the Reynolds numbers and wave numbers that were investigated 
(up to Re = 2000). For the turbulent simulations, the boundary 
condition (Eq. 1) is most simply enforced by solving for the dif- 
ference between the flow and the laminar flow U\\ this difference 
satisfies homogeneous boundary conditions at y — 0 and y — * oo. 

The periods are A, = 0.12L and A t = 0.06 L (about 606; and 
306;). The resolution is as follows. In the y-direction there are 
10 collocation points within about 10 wall units of the wall. In 
the z- and z- directions the spacing is about 20 and 7 wall units, 
respectively. Random numbers are used as initial disturbances 
in order to break any symmetries in the flow field. Their am- 
plitude has to be large enough; otherwise the flow remains lam- 
inar (recall that it is linearly stable). The initial transient takes 
about one cycle and is not included in the phase averages. The 
parity under the transformation <t> — ► <j> + ir is used to double 
the sample (averages are of course taken over z and z as well). 



Fig. 1. Peak wall shear stress vs. Reynolds number. — laminar 
flow; - - - Eq. (14); • direct simulations; o turbulence model. 


Figure 1 shows r m4t / P* as a function of Re. The laminar so- 
lution and the curves corresponding to several values of B in Eq. 
(14) are shown. Dirod simulations were conducted at 
Re - 600, 800, 1000, and 1200. The trend from 800 to 1200 
is steeper than predicted by the theory. The range of Reynolds 
number is too narrow for one to conclude whether the theory is 
invalidated by the results, or whether the disagreement is simply 
a low- Reynolds- number effect. The results obtained with the tur- 
bulence model are included in the figure, but will be discussed 
later. 




Fig. 2. Freestream velocity and wall shear stress vs. phase 
angle. — Uqo/Uq; •• • r in laminar flow; - - - r in turbulent flow, 
a) Re = 600; b) Re = 800. 

At Re = 600, the peak wall stress is close to the laminar value, 
but the flow is not laminar. Figure 2(a) compares the turbulent 
wall stress to the laminar solution; they differ significantly dur- 
ing the deceleration phase. Reversal occurs much later in the 
turbulent flow. The average power consumption per unit area, 
normalized by P<J , is about 0.98 x 10“ 3 in the turbulent flow com- 
pared with 0.83 x 10~ 3 in the laminar flow. The flow is clearly 
not laminar, but is also far from being well -developed turbulence. 
The velocity profile never displays a log layer and the anisotropy 
of the Reynolds- stress tensor is abnormal. Monkewitz and Bun- 
ster (8) visually observed finite- amplitude disturbances, also near 
the phase of reversal, starting between Re ~ 554 and 647. 

Between Reynolds number 600 and 800 a second transition oc- 
curs. Figure 2(b) shows that tbe wall stress is now high dur- 
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ing the high -freest ream -velocity phase. It now has a sharp peak 
with a value much larger than in the laminar flow, and the power 
comsumption is 86% larger than in laminar flow. Strong tur- 
bulence develops around <p = 2.8; a log layer forms in the ve- 
locity profile near <t> = 0 and is maintained well into the de- 
celeration phase, around 4> = 0.7. Reversal now occurs near 
4> = 1.25. The results at Re = 1000 are similar, but the peak in 
the wall stress is much wider. Turbulence during the peak-velocity 
phase was observed in experiments by Merkli and Thomann (9) 
but their Reynolds number was much lower, about 300. See 
also Hino, Sawamoto and Takasu’s experiments in a pipe (10). 



Fig. 3. Wall shear stress vs. phase angle. Re = 1000. 
T / Ttnati - - - cos{<p)\cos{<t>)\. 



Fig. 4. Mean velocity vs. distance from the wall. Re = 1000. 
4 = 0;---* = 2.44; • • • log[j/ + ]/0.41 + 5. 

The agreement with Eq. (12) is tested in Fig. 3. One expects 
T /r mat to have the same shape as cos(*)| cos(*)|, only shifted to 
the left. This is observed only for <p between about -0.3 and 1. This 
is also the range of phases in which a log layer is observed. Figure 
4 shows a log layer at <p = 0, but not at * = 2.44. The profile at 
4> = 2.44 is laminar-like: U + [y+] is far above the usual log layer. 
During the acceleration phase, the favorable pressure gradient de- 
lays the development of the turbulence considerably. Transition 
does not occur until the pressure gradient falls below about 0.02, in 
wall units based on the current wall shear stress, even though the 
near-wall layer is subjected to strong disturbances from the tur- 
bulence remaining farther from the wall. Thus one can attribute 
the lack of agreement with Eq. (12) during that phase to low- 
Revnolds- number effects. This lack of agreement and, the narrow 


range of Reynolds numbers make an accurate check of Eq. (11) im- 
pOS5>ble; one can just note that the value of A is of the order of 3. 



Fig. 5. Profiles in Re = 1000 flow, * = 0. a) Velocity U JU^\ — 

turbulent, laminar, b) Shear stress r/Uj ; total stress, 

- - - Reynolds stress, c) turbulence intensities u' /Ua, etc.: — u', 
v 1 , • * ■ w' . 

Other aspects of the flow will now be described using the 
Re = 1000 case, for which the best sample (four half- cycles) hap- 
pens to be available. Figure 5 shows at * = 0 the mean velocity, 
the shear stress, and the turbulence intensity in each component. 
The near-wall behavior is very similar to that of the usual bound- 
ary layers, but away from the wall the shear stress crosses 0 and 
returns to it only slowly, and the intensities have very long “tails”. 
They extend far heyond the region with significant mean shear U y . 
This is caused bv the lack of entepinment. Notice also how a local 
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minimum in the turbulence intensity, a zero crossing of the shear 
stress, and a zero crossing of the mean shear l\ all coincide near 
y;b> =8. As a rule, the zero crossings of the stress and the shear 
tend to follow each other. Since they imply that the turbulent- 
energy production is 0, they explain the minimum in the energy. 


S-j 






V 





<p 


Fig. 6. Zero crossings of U - U 0 0 vs. phase angle. Re = 1000. 



Fig. 7. Peak value of various quantities over a cycle, vs. y. 

Re = 1000. —(V- Uvo); 100 x r; - 10 x (u' 2 + v 2 + w 2 ). 

Normalized by U 0 . 


The results further show that the structure of the turbulence 
(for instance the anisotropy of the Reynolds-stress tensor) at dif- 
ferent values of y/6 is quite similar. The zero crossings of the stress 
r > of the mean shear U yy and of the deviation from freestream 
U - Uoo all propagate away from the wall at roughly constant ve- 
locity, except close to the wall. Figure 6 shows that the U = 
point starts from the wall at <t> - tt/ 2 with a slope of 1 (the 
laminar value), and that near <f> = 2.5 its velocity abruptly in- 
creases to about 6. The peak values of r, the turbulent energy, 
and U - Uoo over a cycle also decay exponentially as functions 
of yfb (Fig. 7). The decay rate is about 7 times smaller than 
in laminar flow. Thus there is a strong analogy with the laminar 


flow (Eq. 2); the “eddy viscosity” would be uniform at about 40 
times the molecular viscosity (at Re — 1000). 


These findings were somewhat surprising. It was expected that 
the propagation velocity of the waves would scale with the local 
velocity scale, and therefore decrease rapidly away from the wall. 
For the peak values of U - U «, and r, an algebraic decay seemed 
more likely. Note that a constant viscosity scale and an expo- 
nential decay of the velocity scale imply an exponential growth 
of the length scale. On the other hand, it was found that in the 
Reynolds-stress budgets, the role of the various terms did change 
with y: near the wall the usual terms dominate but away from 
it, the turbulent-diffusion term becomes much stronger. Thus the 
similarity is not complete. Even so, the behavior in Figs. 6 and 7 
ought to be explained. 

To further document the behavior of the turbulence the k — e 
model was applied to the flow. The high- Reynolds number version 
of the model was used with standard values for the constants 
(C M = 0.09, C tl = 1.45, C t2 = T9, <r k = 1, <r t = 1.3) (U). 
Only the region y > 26/ was considered: boundary conditions for 
U , k and t at y = 26/ were taken from the direct simulations ( Re = 
1000). Thus near- wall corrections were not needed. For large y 
trivial conditions were used(U - U fc = e = G). The fc-c results 
were in very good agreement with the direct simulations, including 
the uniformity of the eddy viscosity. The agreement was further 
improved by a slight adjustment of C ft1 to 0.085. The direct- 
simulation results also suggested lower values for <7*, around 0.8, 
but such values did not alter the results noticeably. This exercise 
shows that, excluding the wall region, the oscillating boundary 
layer is within the reach of the high-Reynoids-number k — t model. 
Thus, near-wall corrections to the model could be tested with 
confidence that the outer region is treated properly. 


ALGEBRAIC-TURBULENCE-MODEL RESULTS 


The eddy- viscosity model used in this paper is algebraic in the 
sense that the eddy viscosity in a given profile is a specified func- 
tion. However, one of the parameters depends on the phase angle 
(in other problems it could depend on time or on the streamwise 
coordinate) and is determined from solution of an ordinary differ- 
ential equation. The development of the model follows the proce- 
dure originated by Johnson and King (12)i but differs in detail. 
The structure parameter is used to convert the turbulent -energy 
equation to an equation for r, which is integrated over the profile 
to obtain an ordinary differential equation. Models of this type 
can exhibit transition and relaminarization properties in contrast 
to conventional algebraic models. A complete description will be 
given elsewhere. 


The peak wall stress obtained with the model between Re ~ 800 
and 10,000 was shown in Fig. 1. Where the direct-simulation and 
turbulence- model results overlap, the agreement is acceptable. In 
that range both sets of results disagree with Eq. (14) since they do 
not follow a B = constant curve. When Re exceeds about 4.000, 
the agreement between the model and Eq. (14) is seen to improve. 
This suggests that the discrepancy for Re of the order of 1000 is 
indeed due to low- Reynolds- number effects, for which the theory 
is not expected to account. The value of B would then be between 
0.3 and 0.5. On the other hand the Reynolds-number-similarity 
hypothesis and the law- of- the- wall hypothesis, which are central 
to the theory, are also “built into” the turbulence model. In that 
sense, the model does not provide a truly independent check of 
the theory. It is still very useful in assessing the low-Reynolds- 
number effects, which the direct simulations will be unable to do 
in the near future because of their excessive cost. Note also that 
the model predicted relaminarization, between Re = 800 and 600. 
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Fig. 8. Wall shear stress vs. phase angle. Re = 1000. direct 

simulation; turbulence model. 

Figure 8 compares the evolution of the wall stress at 
Re = 1000, computed by direct simulation and with the model. 
The agreement is very good; the model even predicts the quiescent 
phase around 0 — 2 and the surge after <p = 2.5 (although the surge 
occurs slightly earlier with the model). The agreement away from 
the wall (velocity profiles, etc.) is also good. Finally, Fig. 9 shows 
T/ T m*z versus 0 at Re - 10,000. The curve follows Eq. (12) well 
for a much longer part of the cycle than at Re = 1000 (roughly, 
from 0=-lto0 = l). The phase shift 0 O is also smaller. 
The predictions of the model are quite consistent with the theory. 



Fig. 9. Wall shear stress vs. phase angle. Re = 10,000. 

T ’/ T mor i * — cos (0)1 cos( 0 ) | . 

CONCLUSIONS 

The oscillating boundary layer displays a complex behavior, 
both as a function of phase angle and of Reynolds number. Al- 
though it is believed to be linearly stable, it exhibits a first tran- 
sition to a “pre-turbulent” state just below a Reynolds number 
of 600. A second transition, between 600 and 800, allows it to 
generate well-developed turbulence during at least part of the cy- 
cle. During that part of the cycle, it contains a log layer and 
agrees with other aspects of a high- Reynolds- number theory that 
waa also presented. Thus the theory is confirmed to some extent 
by the direct -simulation results. These results even suggest that 
the theory could be simplified further, which does not seem to be 
justified by traditional arguments but is also supported by pre- 
dictions of the k - t model. A new algebraic turbulence model, 
designed and calibrated on the present flow, yielded satisfactory 
agreement with the direct simulations and with the theory in spite 
of the complexity of the flow. 
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